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(SI 8#J8D CD^^JR^O-C. B-X^5-f>B8^B 
4l1 (x) £ TB, (x) J ±BS!2-fS. . * 
n-I 

F(x) = .^ PjBi(x) 

[0 0 2 2 ] ±5£ (3) O^p^^SCitCi 

[0023] -ect, #l&p,tco^-c©jt5i:*fis;^ 

n-l 

F(Xj) = .2^ P i B i (x j >=y j , j=0, •• 

[0 0 2 5 ] *^p,CDSt«n + 3fS>6©tC*J-L. i 
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[0 06 9 ] t^-f>^'J >^(cj:^r^c3n^ 

^iOfetS^Ojg^ttfeOitt* (^9#^) . 
[0 07 0 ] 7 -f ^U* 'J>^L/fcf-^^6, fflreKCffl 

&ft (Hi o«JH) o fit, ±IBffi*CDS&B4 1« 
WO/t^f*ffl^t, x^-f>H«*ifiW4 (H 

1 1 #JI) c 

[oo7i]i«K«i*f-i'B, y->^y 

>^i7^^ y>y<0Saffl*JSfiLfcC4(cJ:o"C* 
CD^^^^a bOtZtiXtoK), ^sEsnsx^ 

[0 07 2] L^O, X^9-f>B»tft«EtfclKC» 

y'vyVT-Z) 4; -eofig (if>^y>^^> 
h) «cfe»*^^^^>wa<DiiK:»JS (Rtt) 
(HI 2#SS) 0 

[0 07 3] CCDm^MLX, BM^^' *>^y> 
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* ttH^- *Kl»r&»SS4:*a&*MStHjt*a I 0 4 
fflA«£JiSmi 0 3-c4ffi5^77^i 

Sioei, K3fc««#ai 0 4*cJ:0 5fot>6 

(t 1( x,), i = 0, n£*5<, Ctitt, 

y y > 1 1 tctew & zftimtfi x^i^ci 

[0 07 9] *-^t>^U>y#Rl 0 ltt. #y 

20 W (X7-7^20 1) B t&^ k£g&cLT*<h 
y v Y t = kdfc*«-&fiia k £, t i St<t 1# i 
fr£ j fc*fl/T:fc5£ (10) T**5. 
[0 08 0] 
[R8] 



t - t 



J 



* j+1 

d (n d 



~ t 



(10) 



[0 0 8 1 ] m.W(DU&±, 0 £ t . , t n ^n 

l*jEa«tt) ^*3<o *fc. t <t.r*st = kdOc 
»ltBa k =x 0 . t ^ t„r*St = kdtc*fbt:« 

[0 0 8 2] 7Y;V^y>^fSl 0 2&, 
> ^i;>^gl o 1 T^fc/cS^* hkd&C>tf-f&ffi 
{a k } GC*tL/T\ *T $ 0 < te? a )l5W&t com 

mmn (*w&a«i#> tfTS? 1 ^*^ * ^* y>^ 

1 2 1 

b v =— 7r~ a ^-i + — 



6 



3 



e 



l k+l 



(11) Ccj^TROCMBbfciilWTS. 
(1 1) fcfc^TA^-^Ja^i. a k , a k + 1 te#£ 
#S*i4«Rl/6 f 2/3, l/6**:7-f^*«R* 

[0 08 3 ] 

[Rg] 

-Ui> 



[0084]^, 7^;^y>^sio2tfe»n 40 

[0 0 8 5] & 0 b k om&¥&>J£it>tcffiT£<DfB. 

v hcr>fflPSd*d/2-Cg#jft^-S ?:/2 0 

4) o 

[0086] — o^ctwna. tnoiaiH£ 



[0 0 8 7] -rttte-fc. ^^2 0 5{C*J^t, 
iSM^©10 3(l y-vb (kd) 6c»-TS7-r;l/ 

CD^, X^^>Eft8lF(t:)£, (kd) 4r-/^ h i L 

(i3), (i4), (is) zm< 0 

[0 08 8] 
[RIO] 



(8) 
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F<t)=2L?p k B k <t) 
nd-1 

F(kd)=f S P| B j(kd)-b k , k = 0, -. n d 
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d 2 F 



d x' 



d 2 F 



nd-l d^B 



t=0 



acrl 



d x 4 



d*B; 



= 0 .-(14) 



t=0 



t=n d d *-* *dx J 



t=n d d 



= 0 -(15) 



[0 0 8 9] mOSv hOP^Pid«:^fS-rSB-X7'7 

(P,) B. *&flnJltLTi»<. 
y 9 r R9RI d -> fc*§£ (^f^204 ©*SSI) 
{C«, m<DWfr\ (Array) SCteffi-T*,, 
[0 090] SI^Sfg#f8 10 411. &]&LtcX A 
>ft«F(t)©£-t,{c:fcW4ffiFct,) tx.iOid 
x.SrtriTf* (^fv^2 06) . 
[0 09 1] ZLX> 11 dx,©**;**!^* (X^ 
5r7'2 0 7) „ 

[0 09 2] XTv72 0 7<DmW&m<,Ci5^X> Id 
x,^?iit)Saf>fc^<D(il (8K£iI) VLTV&tili, A 
^7=-$*sfflrat?^/cCi*>6. ta^^Sl 0 6tCJ:& 

t-^209) . -e^-C&Wft&S. glMdx.SrSr/cft 
f-^x.iLt, twt- 9->^y>^#«l 0 ltC 

(^f -^208) . 
[0 093] HJ^^S 1 0 6 (J, flNltttt&S^Bt 1 0 
3fC<i:D*s£34lfc. d=l. 1/2. l/2\ 

ttiTj-rs (xf-^z 0 9). 

[0 094] 

[lores 2] *sigom2osiss^«s*WTK:siHj-r 

[0 095] l33^#ffSLT:. *jUSJg«gtCfcl»Ttt. 
A#3M9:3 0 5IJ1 ^Tt^rartcD.^t , K*f-r5"<^ r 
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CLAIMS 



[Claim(s)] 

[Claim 1] By repeating the process which generates the approximation curve of the numeric data by 
which the input was carried out [ aforementioned ], in case the spline curve for being the method of 
generating a spline curve from the inputted numeric data, and interpolating numeric data is generated 
While controlling so that the error of the numeric data and the aforementioned approximation curve by 
which the input was carried out [ aforementioned ] approaches 0 In case the aforementioned 
approximation curve is generated, over sampling technique of the aforementioned error of the numeric 
data by which the input was carried out [ aforementioned ] or the numeric data by which the input was 
carried out [ aforementioned ], and the aforementioned approximation curve is carried out. this value by 
which over sampling technique was carried out — filtering processing — carrying out — new data — 
generating — this — the generation method of the spline curve characterized by what is asked for the 
spline curve which interpolates new data 

[Claim 2] The generation method of the spline curve according to claim 1 to which the numeric data by 
which the input was carried out [ aforementioned ] is characterized by the bird clapper from what was 
expressed as a scalar value over the point in 1 -dimensional space. 

[Claim 3] The generation method of the spline curve according to claim 1 to which the numeric data by 
which the input was carried out [ aforementioned ] is characterized by the bird clapper from what was 
expressed as a vector value over the point in I -dimensional space. 

[Claim 4] The generation method of a spline curve given in any 1 of the claims 1-3 characterized by 
interpolating the point of two ****** of the numeric data by which the input was carried out 
[ aforementioned ] with a predetermined primary function, and performing the aforementioned over 
sampling technique. 

[Claim 5] The generation method of a spline curve given in any 1 of the claims 1-3 characterized by 
using the value of B-spline function in a knot for the aforementioned filtering processing. 
[Claim 6] By repeating the process which generates the approximation curved surface of the numeric 
data by which the input was carried out [ aforementioned ], in case the spline surface for being the 
method of generating a spline surface from the inputted numeric data, and interpolating numeric data is 
generated While controlling so that the error of the numeric data and the aforementioned approximation 
curved surface by which the input was carried out [ aforementioned ] approaches 0 In case the 
aforementioned approximation curved surface is generated, over sampling technique of the 
aforementioned error of the numeric data by which the input was carried out [ aforementioned ] or the 
numeric data by which the input was carried out [ aforementioned ], and the aforementioned 
approximation curved surface is carried out. this value by which over sampling technique was carried 
out — filtering processing — carrying out — new data — generating — this — the generation method of the 
spline surface characterized by what is asked for the spline surface which interpolates new data 
[Claim 7] The generation method of the spline surface according to claim 6 to which the numeric data 
by which the input was carried out [ aforementioned ] is characterized by the bird clapper from what was 
expressed as a scalar value over the point in two-dimensional space. 
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[Claim 8] The generation method of the spline surface according to claim 6 to which the numeric data 
by which the input was carried out [ aforementioned ] is characterized by the bird clapper from what was 
expressed as a vector value over the point in two-dimensional space. 

[Claim 9] The generation method of a spline surface given in any 1 of the claims 6-8 characterized by 
interpolating three ****** of the numeric data by which the input was carried out [ aforementioned ] 
with a predetermined primary function, and performing the aforementioned over sampling technique. 
[Claim 10] The generation method of a spline surface given in any 1 of the claims 6-8 characterized by 
using the value of B-spline function in a knot for filtering processing. 

[Claim 1 1] Generation equipment of a spline curve characterized by providing the following. The over 
sampling technique means which carries out over sampling technique of the input data inputted through 
the input means, and carries out the generation output of the new data. A filtering means to input the 
data outputted from the aforementioned over sampling technique means, and to perform predetermined 
filtering processing to these data. A interpolation curvilinear generation means to generate the spline 
curve for interpolating the data outputted from the aforementioned filtering means. An error calculation 
means to calculate the error over the aforementioned input data of the spline curve generated with the 
aforementioned interpolation curvilinear generation means, and an output means to change and output 
the spline curve group generated with the aforementioned interpolation curvilinear generation means to 
a single spline curve. 

[Claim 12] Generation equipment of a spline curve according to claim 1 1 with which the 
aforementioned input data inputted through the aforementioned input means is characterized by the bird 
clapper from the scalar value over the point in 1 -dimensional space. 
[Claim 13] Generation equipment of a spline curve according to claim 1 1 with which the 
aforementioned input data inputted through the aforementioned input means is characterized by the bird 
clapper from the vector value over the point in 1 -dimensional space. 

[Claim 14] Generation equipment of a spline curve given in any 1 of the claims 1 1-13 to which the 
aforementioned over sampling technique means is characterized by interpolating the point of two 
****** of the aforementioned input data with a predetermined primary function, and performing over 
sampling technique. 

[Claim 15] Generation equipment of a spline curve given in any 1 of the claims 11-13 characterized by 
for the aforementioned filtering means making the value of B-spline function in a knot a filter factor, 
and carrying out filtering. 

[Claim 16] Generation equipment of a spline surface characterized by providing the following. The over 
sampling technique means which carries out over sampling technique of the input data inputted through 
the input means, and carries out the generation output of the new data. A filtering means to input the 
data outputted from the aforementioned over sampling technique means, and to perform predetermined 
filtering processing to these data. A interpolation curved-surface generation means to generate the spline 
surface for interpolating the data outputted from the aforementioned filtering means. An error 
calculation means to calculate the error over the aforementioned input data of the spline curve generated 
with the aforementioned interpolation curved-surface generation means, and an output means to change 
and output the spline-surface group generated with the aforementioned interpolation curved-surface 
generation means to a single spline surface. 

[Claim 17] Generation equipment of a spline surface according to claim 16 with which the 
aforementioned input data inputted through the aforementioned input means is characterized by the bird 
clapper from the scalar value over the point in two-dimensional space. 
[Claim 18] Generation equipment of a spline surface according to claim 16 with which the 
aforementioned input data inputted through the aforementioned input means is characterized by the bird 
clapper from the vector value over the point in two-dimensional space. 

[Claim 19] Generation equipment of a spline surface given in any 1 of the claims 16-18 to which the 
aforementioned over sampling technique means is characterized by interpolating three ****** with a 
predetermined primary function mutually [ the aforementioned input data inputted through the 
aforementioned input means ], and performing over sampling technique. 
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[Claim 20] Generation equipment of a spline surface given in any 1 of the claims 16-18 characterized by 
for the aforementioned filtering means making the value of B-spline function in a knot a fill coefficient, 
and carrying out filtering. 

[Claim 21] Generation equipment of a spline curve'given in any 1 of the claims 11-13 characterized by 
inputting the numeric data which the aforementioned input means becomes from the digital signal which 
samples an analog signal and is obtained as the aforementioned input data. 

[Claim 22] Generation equipment of a spline surface given in any 1 of the claims 16-18 characterized by 
inputting the numeric data which the aforementioned input means becomes from the digital signal which 
samples an analog signal and is obtained as the aforementioned input data. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[The technical field to which invention belongs] this invention relates to the generation method and 
generation equipment of the spline curve from the numeric data obtained by especially measurement, 
and a curved surface about the generation method of the spline curve and curved surface which are used 
in a computer graphics or CAD (computer support die design). 
[0002] 

[Description of the Prior Art] An analog signal is measured and the Prior art which changes into a spline 
curve or a curved surface the numeric data which changes this into a digital signal etc. and is obtained is 
explained with reference to drawing 21 from drawing 13 . 

[0003] First, from drawing 13 to drawing 16 is explained with reference to an outline. 
[0004] In order to presume the value 1302 in the position (point which is not sampled) where the value 
is not given to the numeric data 1301, such as observed value, such as atmospheric temperature and 
atmospheric pressure, and a digital signal which sampled the analog signal, or to see the distribution of 
the whole value and to graph-ize, the need of asking for the function which interpolates numeric data 
often starts. 

[0005] Moreover, two-dimensional [ objective ] or a 3-dimensional image is measured using a 
measuring instrument, the obtained sequence of points are changed into a curve or a curved surface, and 
there is also a demand used for the design of a product model. 

[0006] As the interpolation method which fills this demand most simply, as shown in drawing 14 , there 
is a method of connecting the sampling data (numeric data) of two ****** in a straight line. However, 
by this method, the configuration of the function 1401 obtained does not become smooth. Drawing 14 
shows 1401, as a result of connecting the numeric data 1301 of drawing 13 in a straight line. 
[0007] Since it is thought that the inputted numeric data was sampled from the smooth signal (for 
example, continuous analog signal) in many cases, the function (1401 of drawing 14 ) generated by 
linear interpolation will become an original function (function showing a former signal) and the original 
thing which made the mistake in becoming. 

[0008] On the other hand, the method of interpolating using a polynomial is one method for connecting 
sampling data smoothly. Drawing 15 shows the curve which connected the numeric data 1301 of 
drawing 13 by the polynomial. However, the order of a polynomial becomes high, and since the 
polynomial 1501 generated has an intense vibration, it is not suitable for practical use by this method, in 
many cases, as the number of data increases. 

[0009] Then, the spline function which is the partition-polynomial which connected a different 
polynomial for every section is used well. 

[0010] In a part for the connection of the polynomial of a spline function, it has a continuity according 
to the order of a polynomial to be used. For example, when using the 3rd polynomial, in a part for a 
connection, the derivative to the second floor has a continuity. For this reason, since the conditions 
about a continuity are loose when it interpolates by one polynomial, the smooth function 1601 with little 
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(1) 



vibration is obtained. 

[001 1] Although what is necessary is just to decide the polynomial of each section in order to perform 
interpolation using the spline function, asking for the coefficient of a polynomial directly is not desirable 
in respect of a calculation error. 

[0012] Then, usually, a spline function is expressed as linear combination of B-spline function, and the 
simultaneous equations about the coefficient are solved. 

[0013] Since B-spline function has a local base, it is restricted near the diagonal line that the element 
which is not 0 in the coefficient matrix of simultaneous equations appears, and a solution is calculated 
stably and at high speed. 

[0014] Next, the concrete calculation method is explained with reference to drawing 21 from drawing 
17. 

[0015] The inputted numeric data is set with (xi, yi), xi<xi +1, i= 0, — , n by making n into a positive 
integer (refer to drawing 17 ). Although the connection of the polynomial of a spline function is called a 
"knot" or "joint", below, it has each xi as a knot and how to ask for the spline fiuiction which 
interpolates the inputted numeric data is explained. Moreover, only a cubic-spline function is explained 
on account of explanation. 

[0016] First, a knot is extended to right and left. That is, it sets with x-3=x-2=x-l=x0 and 
xn+3=xn+2 : =xx+l= : xn. At this time, the following (m-1) B-spline functions Bm and j (x) are defined by 
the following formula (1) and (2) to a knot vector (x-3, x-2, --, xn+3). 
[0017] 
[Equation 1] 

B t j(x)=l (Xj2Sx<x j + 1 ) , 0 (*<Dtf!) 

B ^ (i>= 6^ B -- i - (x)+ 4S. B »- i - j+i<,i> - (m " ) " <2> 

[001 8] However, when 0 appears in a denominator, the term is set with 0. 

[0019] Here, since only the 3rd spline function 1801 (refer to drawing 18 ) is treated, the B-spline 
functions B4 and j (x) are written as "Bj (x)." 

[0020] Then, spline-function F (x) is expressed like the following formula (3) as linear combination of 

the B-spline function Bj (x). 

[0021] 

[Equation 2] 
n-i 

F(x)=.22 p |Bj(x) «.(3) 

[0022] By asking for the coefficient pi of an upper formula (3), the spline function for which it asks 
becomes settled. 

[0023] Then, the simultaneous equations about a coefficient pi are stood. First, in each knot, if the value 

of input data and the value of a spline function are equal, it will set (refer to the following formula (4)). 

[0024] 

[Equation 3] 
n-1 

F(Xj) = .^: PjBi(Xj)=yj, j=0, — , n -(4) 

[0025] To the number of Unknowns pi being n+3, since the number of input data is n+1, it runs short of 
two conditions by the upper formula (4). Then, the following formula (5) which are the endpoint 
conditions that the derivative of the second floor in ends is 0, and (6) are added. 
[0026] 
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d x 



n-1 



d 2 B : 



d x 



= 0 



(5) 



x 0 



d^F 



d x' 



n-1 



d 2 B 



d x 



= 0 



(6) 



x=x r 



[0027] pi can be found if an upper formula (4), (5), and (6) are solved as simultaneous equations. 

[0028] Also about the case more than two-dimensional, it can come back and solve to one dimension. 

Since the same is more said of high order origin, a two-dimensional case is explained. 

[0029] The inputted data shall be located in a line in the shape of a grid. That is, it considers as xi, i= 0, - 

-, m and yj, j= 0, — , the thing by which the value zij is given to (xi, yj) to n (refer to drawing 19 ). 

[0030] If extension of a knot is performed like [ direction / each / x directions and / of y ] the case of one 

dimension and the 3rd spline function is expressed using B-spline function, it will be given by the 

following formula (7). 

[0031] 

[Equation 5] 

■-I n-1 



F (x, y) = 



Sj5 



P ij B i 



(x) Cj <y) 



(7) 



[0032] However, Cj (y) is B-spline function corresponding to yj. 

[0033] At this time, interpolation in one dimension is first performed to each i. That is, from input data 
zij, the function Fi of Variable y (y) is carried out like [ in the above-mentioned one dimension ], and it 
asks for it (refer to drawing 20 ). 
[0034] 
[Equation 6] 



n-1 



i 



(y) =: 



C= (y) 



<8) 



[0035] It performs 1 -dimensional interpolation each about j shortly to the solution qij of an upper 
formula (8). That is, it asks for the function Fj of Variable x (x) by the method in one dimension from qij 
(refer to drawing 21 ). 
[0036] 
[Equation 7] 



[0037] pij obtained by the upper formula (9) serves as a coefficient of B-spline-function expression of 

the two-dimensional spline function for which it asks. 

[0038] 

[Problem(s) to be Solved by the Invention] The trouble of the interpolation by the above-mentioned 
conventional spline function is explained below with reference to drawing 22 and drawing 23 . 
[0039] By the interpolation method by the above-mentioned conventional spline (B-spline) function, 
vibration of the function generated is suppressed as compared with interpolation by the polynomial. 
[0040] However, it may still change a lot to a slight change of input data. 

[0041] This phenomenon is explained about the case where the position of data is unequal. For example, 
if all the values of input data are the same values as shown in drawing 22 , naturally the result of 
interpolation will become the straight-line- like function 2201. 
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[0042] However, when change is included in the section where the function of data is short, by the 
spline function 2301 (refer to drawing 23 ) generated, the error is amplified, and if the value of a 
function 2201 (refer to drawing 22 ) and a function 2301 is compared in the section which the interval of 
data is opening, it will come out making a mistake in being big (error). 

[0043] Moreover, it is decided by the inputted numeric data, and the position of the knot of a spline 
function will be made at equal intervals, when a knot at equal intervals is desirable. 
[0044] Moreover, a spline function cannot be decided to be a meaning unless two-dimensional data are 
distributed in the shape of a grid. 

[0045] therefore - while also being able to suppress change of the function out of which succeeds in this 
invention in view of the trouble of the above-mentioned conventional technology, and change of input 
data comes only and which is generated in a certain case in few amount the knot of a spline function was 
inputted - it aims at offering the generation method and equipment of the spline curve or curved surface 
which makes it possible to be able to choose freely, without being dependent on numeric data 
[0046] 

[Means for Solving the Problem] In order to attain the aforementioned purpose, this invention is the 
method of generating a spline curve from the numeric data which measured the analog signal and was 
obtained. It carries out by repeating the step which generates the approximation curve of the numeric 
data into which generation of the spline curve which interpolates numeric data was inputted, and 
bringing the error from numeric data close to 0. At the step which generates an approximation curve, 
over sampling technique of the aforementioned numeric data or the aforementioned error is carried out, 
it filters further, new data are generated, and it asks for the spline curve which interpolates this generated 
data. 

[0047] What was expressed as a scalar value over the point in 1 -dimensional space as the 
aforementioned numeric data can be taken. 

[0048] Or what was expressed as a vector value over the point in 1 -dimensional space as the 
aforementioned numeric data can also be taken. 

[0049] In the generation method of the spline curve of this invention, over sampling technique can be 
performed by interpolating two adjoining points with a primary function. 

[0050] Furthermore, the value of B-spline function in a knot can be used for filtering in the generation 
method of the spline curve of this invention. 

[0051] The generation method of the spline surface of this invention is the method of generating a spline 
surface from the numeric data which measured the analog signal and was obtained. It carries out by 
repeating the step which generates an approximation curved surface to the numeric data into which 
generation of the spline surface which interpolates numeric data was inputted, and bringing the error 
from numeric data close to 0. It asks for the spline surface which carries out over sampling technique of 
the aforementioned numeric data or the aforementioned error, filters further, generates new data, and 
interpolates the generated data at the step which generates an approximation curved surface. 
[0052] What was expressed as a scalar value over the point in two-dimensional space as the 
aforementioned numeric data can be taken. 

[0053] Or what was expressed as a vector value over the point in two-dimensional space as the 
aforementioned numeric data can also be taken. 

[0054] In the generation method of the spline surface of this invention, over sampling technique can be 
performed by interpolating the value over three adjoining points with a primary function. 
[0055] Furthermore, the value of B-spline function in a knot can be used for filtering in the generation 
method of the spline surface of this invention. 

[0056] An input means by which the generation equipment of the spline curve of this invention inputs 
the numeric data which measured the analog signal, An over sampling technique means to carry out over 
sampling technique of the numeric data inputted with the aforementioned input means, and to generate 
new data, A filtering means to filter the data generated with the aforementioned over sampling technique 
means, A interpolation curvilinear generation means to generate the spline curve which interpolates the 
data filtered with the aforementioned filtering means, It consists of an error calculation means to 
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calculate the error over the aforementioned numeric data of the spline curve generated with the 
aforementioned interpolation curvilinear generation means, and an output means to change and output 
the spline curve group generated with the aforementioned interpolation curvilinear generation means to 
a single spline curve. 

[0057] The aforementioned input means can input the scalar value over the point in 1 -dimensional space 
as numeric data. 

[0058] Or the aforementioned input means can also input the vector value over the point in 1- 
dimensional space as numeric data. 

[0059] In the generation equipment of the spline curve of this invention, over sampling technique can be 
performed by interpolating the value over two points which the aforementioned over sampling technique 
means adjoins with a primary function. 

[0060] Furthermore, in the generation equipment of the spline curve of this invention, the 
aforementioned filtering means can use the value of B-spline function in a knot for filtering. 
[0061] An input means by which the generation equipment of the spline surface of this invention inputs 
the numeric data which measured the analog signal, An over sampling technique means to carry out over 
sampling technique of the numeric data inputted with the aforementioned input means, and to generate 
new data, A filtering means to filter the data generated with the aforementioned over sampling technique 
means, A interpolation curved-surface generation means to generate the spline surface which 
interpolates the data filtered with the aforementioned filtering means, It consists of an error calculation 
means to calculate the error over the aforementioned numeric data of the spline surface generated with 
the aforementioned interpolation curved-surface generation means, and an output means to change and 
output the spline-surface group generated with the aforementioned interpolation curved-surface 
generation means to a single spline surface. 

[0062] The aforementioned input means can input the scalar value over the point in two-dimensional 
space as numeric data. 

[0063] Or the aforementioned input means can also input the vector value over the point in two- 
dimensional space as numeric data. 

[0064] In the generation equipment of the spline surface of this invention, over sampling technique can 
be performed by interpolating the value over three points which the aforementioned over sampling 
technique means adjoins with a primary function. 

[0065] Furthermore, in the generation equipment of the spline surface of this invention, the 
aforementioned filtering means can use the value of B-spline function in a knot for filtering. 
[0066] 

[Function] The principle of this invention and an operation are explained in detail below with reference 
to drawing 12 from drawing 7 . 

[0067] First, the data which carried out over sampling technique from the inputted data (refer to drawing 
7 ) are generated (refer to drawing 8 ). The interval of the sampling point (sample takeoff point) of the 
data (sampling data) which over sampling technique generally inputted is subdivided, and although it is 
making the data in the subdivided point, it shall be called over sampling technique not only including 
when the sampling period of the data inputted here is equal, but making the data in arrangement equal 
when it is arrangement with an unequal sampling period. 

[0068] Although there are various methods in over sampling technique, after interpolating between data 
with a primary function here, how to calculate the value of the data in a position to ask is explained as an 
example. 

[0069] The data generated by over sampling technique are filtered. Filtering processing obtains new data 
from the original data in quest of the weighted mean about surrounding data. The data generated by 
filtering processing become what has boom hoisting looser than the original data (refer to drawing 9 ). 
[0070] The value in the point corresponding to the knot of the spline function used for interpolation is 
chosen from the filtered data (refer to drawing 10 ). And a spline function is generated using the method 
explained as the above-mentioned Prior art (refer to drawing 1 1 ). 

[0071] When the data used for interpolation went via processing of over sampling technique and 
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filtering, the boom hoisting is made loose and the big change to the spline function generated does not 
appear. 

[0072] However, since the data used when a spline function is generated are over sampling technique 
and data obtained by carrying out filtering processing about input data and are not the input data itself, 
there is a difference (error) in the original input data (sampling data) and the value of the spline function 
in the position (sample takeoff point) (refer to drawing 12 ). 
[0073] It filters with over sampling technique again to this difference. 

[0074] If a processing result appears as a value in a knot at this time, interpolation to an error will be 

performed and the precision of a spline function will be raised. On the other hand, when not appearing 

as a value in a knot, the interval of a knot is reduced and a spline function is generated. 

[0075] Finally interpolation of the request to input data is attained by carrying out this invention in this 

way, and cutting down gradually the error of the data obtained by interpolation, and the original input 

data (sampling data). 

[0076] 

[Embodiments of the Invention] The form of operation of this invention is explained in detail below 

with reference to a drawing. 

[0077] 

[Operation form 1] The 1st operation form of this invention is explained with reference to drawing 1 . 
An input means 105 to input a scalar value [ as opposed to the point in 1 -dimensional space in this 
operation form ] as shown in drawing 1 , The over sampling technique means 101 which carries out over 
sampling technique of the numeric data inputted through the input means 105, The filtering means 102 
which carries out filtering processing of the data which were outputted from the over sampling 
technique means 101, and which were carried out, A interpolation curvilinear generation means 103 
against the new data outputted from the filtering means 102 to generate a spline function as 
interpolation, An error calculation means 104 to search for the error over the numeric data into which 
the spline curve generated by the interpolation curvilinear generation means 103 was inputted, An 
output means 106 to change and output the spline curve group generated with the interpolation 
curvilinear generation means 103 to a single spline curve, The error of the value of the input data called 
for by ****** and the error calculation means 104 and the interpolation data to the sample point is 
supplied to the over sampling technique means 101, and the over sampling technique means 101 carries 
out over sampling technique of this error, and outputs it. 

[0078] Next, processing operation of this operation gestalt is explained with reference to the flow chart 
shown in drawing 1 and drawing 2 . Below, the numeric data of 1 -dimensional one into which n was 
inputted as a predetermined positive integer is set with (ti, xi), i= 0, — , n. This shows that the sampled 
value in a sampling point ti is xi. Moreover, the spline curve generated in the interpolation curvilinear 
generation means 103 sets the interval of the knot with d with a knot at equal intervals. The initial value 
of d is set with 1 on account of explanation here. 

[0079] The over sampling technique means 101 calculates the value in each knot with interpolation by 
the primary function of input data (Step 201). namely, the value [ as opposed to / when making k into an 
integer / knot t=kd ] ak - tj<=t<tj +1 it asks by the following formula (10) from j 
[0080] 
[Equation 8] 

t - t j 

a k = Xj + <x j+1 - Xj ) x - — ...(10) 

t i+l~ t i 

[0081] It sets with 0 <=t0 and tn<=ndd (nd is a positive integer) on account of explanation. Moreover, to 
t=kd which is t<t0, it considers as ak=xn to ak=x0 and t=kd which is t>=tn. 
[0082] The filtering means 102 performs digital filtering processing in which a sum-of-products 
operation (collapsing operation) with a filter factor is performed preferably, to the train {ak} of the value 
over each knot kd for which it asked with the over sampling technique means 101 (Step 202). The value 
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of B-spline function to a knot is used for a filter. Specifically, the new value bk is calculated by the 
following formula (1 1). In the following formula (1 1), input data train ak-1, ak, the coefficients 1/6 by 
which multiplication is carried out to ak+1, 2/3, and 1/6 show the filter factor. 
[0083] 
[Equation 9] 

1 2 1 



b„ =- 



6 



a k-l 4 



3 



+ 



8 



k+1 



(11) 



[0084] Next, the size of the value bk acquired with the filtering means 102 is investigated (Step 203). 
[0085] If the value of bk is below a predetermined value (default value) defined beforehand, since an 
approximate precision does not improve any more, the interval d of a knot will be replaced by d/2 at the 
spline curve whose interval of a knot is d (Step 204). 

[0086] On the other hand, if larger than default value, the interpolation spline curve by the interpolation 
curvilinear generation means 103 will be generated (Step 205). 

[0087] That is, in Step 205, the interpolation curvilinear generation means 103 generates a spline curve 
using the filter result of an operation bk to a knot {kd}. In that case, spline curve F (t) is expressed like 
the following formula (12) as linear combination of the B-spline function Bk (t) which has {kd} as a 
knot, and the simultaneous equations (13) about a coefficient pk, (14), and (15) are solved. 
[0088] 

[Equation 10] 
nd-l 



F(t) = ^ P k B k (t) 
k— 3 * K 



(12) 



nd-l 



F(kd)=?2C P jB j(kd)=b k , k = 0, n d 



-(13) 



d x 



nd-l 



t=0 



dx' 



= 0 



t=0 



(14) 



d 2 F 



dx 



d 2 B 



nd-l u"D| 
t-n d d = ^ Pl d^ 



t=n d d 



= 0 



(15) 



[0089] The coefficient {pi} of B-spline function corresponding to the interval d of the same knot is 
added one by one. On the other hand, when the knot interval d changes (processing of Step 204), it 
stores in another array (Array). 

[0090] every of spline curve F (t) which generated the error calculation means 104 ~ the difference dxi 

of value F (ti) and xi in ti is calculated (Step 206) 

[0091] And the size of Difference dxi is investigated (Step 207). 

[0092] In judgment processing of Step 207, if Difference dxi is below the predetermined value (default 
value) defined beforehand, since the input data was interpolated, a spline curve will be outputted in the 
procedure explained to the following depended on the output means 106 (Step 209). Otherwise, it inputs 
into the over sampling technique means 101 by using Error dxi as the new data xi (Step 208). 
[0093] The output means 106 changes d= 1, 1/2, 1/22 and - which were generated by the interpolation 
curvilinear generation means 103, and the spline curve group corresponding to 1/21. into the spline curve 
whose interval of a knot is 1/21., applies the coefficient pi of B-spline function, and makes and outputs 
one spline curve (Step 209). 
[0094] 



http://www4.ipdl.jpo.go.jp/cgi-bin/tran_web_cgi_ejje 



9/4/2003 



Page 8 of 10 



[Operation gestalt 2] The 2nd operation gestalt of this invention is explained below. Drawing 3 is the 
block diagram showing the composition of the 2nd operation gestalt of this invention. 
[0095] With reference to drawing 3 , the input means 305 inputs the vector value Vi over the point ti in 
1 -dimensional space in this operation gestalt. Each processing of the over sampling technique means 
301, the filtering means 302, the interpolation curvilinear generation means 303, the error calculation 
means 304, and the output means 306 performs the same processing as the over sampling technique 
means 101 in the operation gestalt of the above 1st, the filtering means 102, the interpolation curvilinear 
generation means 103, the error calculation means 104, and the output means 106 about the value 
(namely, scalar value) of each component of a vector value. For this reason, explanation of processing of 
these meanses is omitted. 
[0096] 

[Operation gestalt 3] The 3rd operation gestalt of this invention is explained below. Drawing 4 is the 
block diagram showing the composition of the 3rd operation gestalt of this invention. 
[0097] Drawing 4 is referred to. the 3rd operation gestalt of this invention An input means 405 to input 
the scalar value over the point in two-dimensional space, The over sampling technique means 401 which 
carries out over sampling technique of the numeric data inputted by the input means 405, or the error 
calculated with the error calculation means 404 as input data, A filtering means 402 to filter the data 
calculated with the over sampling technique means 401, A interpolation curved-surface generation 
means 403 to generate the spline function which interpolates the new data calculated with the filtering 
means 402, It has an error calculation means .404 to calculate the error of the spline function and data 
which were generated with the interpolation curved-surface generation means 403, and an output means 
406 to output the data of the generated spline surface. 

[0098] With reference to the flow chart of drawing 4 and drawing 5 , operation of the 4th operation 
gestalt of this invention is explained in detail below. 

[0099] Here, the inputted numeric data shall be located in a line in the shape of a grid (two-dimensional 
matrix) on account of explanation. That is, it considers as ui, i= 0, -, m and vj, j= 0, -, the thing by 
which the value xij is given to the lattice point (ui, vj) to n. 

[0100] And to each grid which consists of four points (ui, vj), (ui+1, vj), (ui, vj+1), and (ui+1, vj+1), it 
shall be underlined with the diagonal line which connects (ui+1, vj) (ui, vj+1), and triangulation shall be 
carried out. 

[0101] Since it can process similarly if triangulation, such as the Delaunay triangulation, is performed 
even when having not stood in a line in the shape of a grid, the detail when having not stood in a line in 
the shape of a grid is omitted. 

[0102] The interval of the knot of a spline function sets the direction of u, and the direction of v with d. 
Initial value of d is set to 1 . 

[0103] The over sampling technique means 401 calculates the value in each knot lattice point with 

interpolation by the primary function of input data (Step 501). That is, when you make k and 1 into an 

integer, let the value ckl in lattice point (u, v) = (kd, Id) of a knot be the value which interpolated the 

value in the triangular vertex where (u, v) are contained with the primary function. 

[0104] Therefore, when setting to (ui, vj), (ui+1, vj), and (ui, vj+1) the vertex of the triangle contained 

first (u, v), it asks for barycentric coordinates bO, bl, and b2 from the following formula (16), (17), and 

(18). 

[0105] bO ui+blui+l+b2 ui=u - (16) bO vj+bl vj+b2vj+l=v - (17) b0+bl+b2=l - (18) [0106] Thus, 
barycentric coordinates bO, bl, and b2 for which it asked It uses, and a value ckl is carried out like the 
following formula (19), and is calculated. 
[0107] 

ckl=b0 xij+bl xi+1, j+b2xi, j+1 - (19) [0108] It sets on account of explanation with 0 <=u0, um<=mdd, 
0 <=v0, and vn<=ndd (md and nd are a positive integer). 

[0109] In the exterior of the grid of input data, if u coordinate is within the limits of a grid, the nearest 
thing in the side of the outermost periphery will be chosen, and u coordinate will consider as the value 
of the same point. Moreover, it treats similarly about v coordinate. 



http://www4.ipdl.jpo.go.jp/cgi-bin/tran_web_cgi_ejje 9/4/2003 



Page 9 of 10 



[0110] If the grid of u coordinate and v coordinate is out of range, let them be a value over the nearest 
thing among the four peaks of the outermost periphery of a grid. 

[0111] The filtering means 402 performs a two-dimensional filter operation (sum-of-products operation 
of a filter factor and an array {ckl}) to the array {ckl} of a value to the lattice point (kd, Id) of each knot 
for which it asked with the over sampling technique means 401 (Step 502). The product of the value of 
B-spline function to the knot lattice point is used for a filter. Specifically, a new value {ekl} is 
calculated by the following formula (20). 
[0112] 

[Equation 11] 

1 1 1 1 4 

e kl~ ^lT c k~l, l-l + ^ c kJ-l+J^ c k+M-l+^ c k-l, l + i" c k, 1 

1111 

+y c k+l, l + lT^ c k-l. l + l+^ c k, l + l + 7f^T c kU. 1+1 - (20) 

[0113] Next, the size of the value ekl acquired with the filtering means 402 is investigated (Step 503). If 
the size of a value ekl is below the predetermined value (default value) defined beforehand, since an 
approximate precision does not improve any more, the interval d of a knot will be replaced by l/2d at 
the spline function whose interval of a knot is d (Step 504). On the other hand, if larger than default 
value, in Step 505, the interpolation spline surface by the interpolation curved-surface generation means 
403 will be generated. 

[0114] The interpolation curved- surface generation means 403 generates a spline function using ekl to 
(kd, Id). If 3rd spline-surface F (u, v) is expressed using B-spline function, it will become the following 
formula (21). 
[0115] 

[Equation 12] 

F (u. v) =jZJ jEZpjjBj (u) Cj (v) —(21) 

[01 16] However, Bi (u) and Cj (v) are B-spline functions corresponding to {kd} and {Id} in a knot, 
respectively. 

[0117] At this time, interpolation in one dimension is first performed to each i. That is, from input data 
eij, the function Fi of Variable v (v) is carried out like [ in the above-mentioned one dimension ], and it 
asks for it. 
[0118] 

[Equation 13] 

Fj (v) =jZ2qjjCj (v) -(22) 

[01 19] It performs 1 -dimensional interpolation each about j shortly to the solution qij of an upper 
formula (22). That is, it asks for the function Fj of Variable u (u) by the method in one dimension from 

qy. 

[0120] 

[Equation 14] 

Bd-1 

Fj (u) =^PijBi (u) -(23) 
[0121] It is the value which pij obtained by the upper formula (23) calculates. 

[0122] The coefficient {pij} of B-spline surface corresponding to the same knot interval d is added one 
by one. When d changes, it stores in another array. 
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[0123] The error calculation means 404 calculates the difference dxij (error) with the values F (ui, vj) 
and xij in ** (ui, vj) of generated spline-surface F (u, v) (Step 506). 

[0124] And the value of dxij is investigated (Step 507). Since the input data was interpolated when dxij 
was below default value, it progresses to the output means 406. Otherwise, it inputs into the over 
sampling technique means 401 by using Error dxij as the new data xij (Step 508). 

[0125] The output means 406 changes d= 1, 1/2, 1/22 and — which were generated with the interpolation 
curvilinear generation means 403, and the spline-surface group corresponding to 1/21. into the spline 
surface whose interval of a knot is 1/21., applies the coefficient pij of B-spline function, and makes and 
outputs one spline surface (Step 509). 
[0126] 

[Operation gestalt 4] The 4th operation gestalt of this invention is explained with reference to drawing 
6 . In the 4th operation gestalt of this invention, the input means 605 inputs the vector value Vij over the 
point in two-dimensional space (ui, vj). 

[0127] Since the same processing as an over sampling technique means 401 to correspond in the 3rd 
operation gestalt about each component of a vector value, respectively, the filtering means 402, the 
interpolation curved-surface generation means 403, the error calculation means 404, and the output 
means 406 is performed, each processing of the over sampling technique means 601, the filtering means 
602, the interpolation curved-surface generation means 603, the error calculation means 604, and the 
output means 606 is omitted for details. 
[0128] 

[Effect of the Invention] It has the effect that change of the function obtained as a result can also 
suppress small to a small change, by having constituted so that the spline interpolation of the data which 
carried out over sampling technique may carry out about the value acquired by carrying out filtering 
processing to the inputted numeric data according to the interpolation method of numeric data and the 
equipment concerning this invention, as having explained above. For this reason, a interpolation 
function smoother than the conventional technique can be obtained. 



[Translation done.] 
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* NOTICES * 

Japan Patent Office is not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



PRIOR ART 



[Description of the Prior Art] An analog signal is measured and the Prior art which changes into a spline 
curve or a curved surface the numeric data which changes this into a digital signal etc. and is obtained is 
explained with reference to drawing 21 from drawing 13 . 

[0003] First, from drawing 13 to drawing 16 is explained with reference to an outline. 
[0004] In order to presume the value 1302 in the position (point which is not sampled) where the value 
is not given to the numeric data 1301, such as observed value, such as atmospheric temperature and 
atmospheric pressure, and a digital signal which sampled the analog signal, or to see the distribution of 
the whole value and to graph-ize, the need of asking for the function which interpolates numeric data 
often starts. 

[0005] Moreover, two-dimensional [ objective ] or a 3-dimensional image is measured using a 
measuring instrument, the obtained sequence of points are changed into a curve or a curved surface, and 
there is also a demand used for the design of a product model. 

[0006] As the interpolation method which fills this demand most simply, as shown in drawing 14 , there 
is a method of connecting the sampling data (numeric data) of two ****** ina straight line. However, 
by this method, the configuration of the function 1401 obtained does not become smooth. Drawing 14 
shows 1401, as a result of connecting the numeric data 1301 of drawing 13 in a straight line. 
[0007] Since it is thought that the inputted numeric data was sampled from the smooth signal (for 
example, continuous analog signal) in many cases, the function (1401 of drawing 14 ) generated by 
linear interpolation will become an original function (function showing a former signal) and the original 
thing which made the mistake in becoming. 

[0008] On the other hand, the method of interpolating using a polynomial is one method for connecting 
sampling data smoothly. Drawing 1 5 shows the curve which connected the numeric data 1301 of 
drawing 13 by the polynomial. However, the order of a polynomial becomes high, and since the 
polynomial 1501 generated has an intense vibration, it is not suitable for practical use by this method, in 
many cases, as the number of data increases. 

[0009] Then, the spline function which is the partition-polynomial which connected a different 
polynomial for every section is used well. 

[0010] In a part for the connection of the polynomial of a spline function, it has a continuity according 
to the order of a polynomial to be used. For example, when using the 3rd polynomial, in a part for a 
connection, the derivative to the second floor has a continuity. For this reason, since the conditions 
about a continuity are loose when it interpolates by one polynomial, the smooth function 1601 with little 
vibration is obtained. 

[001 1] Although what is necessary is just to decide the polynomial of each section in order to perform 
interpolation using the spline function, asking for the coefficient of a polynomial directly is not desirable 
in respect of a calculation error. 

[0012] Then, usually, a spline function is expressed as linear combination of B-spline function, and the 
simultaneous equations about the coefficient are solved. 

[0013] Since B-spline function has a local base, it is restricted near the diagonal line that the element 
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which is not 0 in the coefficient matrix of simultaneous equations appears, and a solution is calculated 
stably and at high speed. 

[0014] Next, the concrete calculation method is explained with reference to drawing 21 from drawing 

17. . 
[0015] The inputted numeric data is set with (xi, yi), xi<xi +1, i= 0, --, n by making n into a positive 
integer (refer to drawing 17 ). Although the connection of the polynomial of a spline function is called a 
"knot" or "joint", below, it has each xi as a knot and how to ask for the spline function which 
interpolates the inputted numeric data is explained. Moreover, only a cubic-spline function is explained 
on account of explanation. 

[0016] First, a knot is extended to right and left. That is, it sets with x-3</SUB> =x-2=x-l=x0 and 
xn+3=xn+2=xx+l=xn. At this time, the following (m-1) B-spline functions Bm and j (x) are defined by 
the following formula (1) and (2) to a knot vector (x-3, x-2, -, xn+3). 
[0017] 
[Equation 1] 

B lfJ (x)=l ( Xj =!Sx<x j + 1 > , 0 (*<DtiL) -(» 

o X ~ Xj B , .( X )+ Xj * a X B m -i i+1 (x).(mS2) -(2) 

B -.J (x) -x j+m _ 1 -x j B «>- 1 .J W * j+1D -x j+1 Bl ' J+1 

[001 8] However, when 0 appears in a denominator, the term is set with 0. 

[0019] Here, since only the 3rd spline function 1801 (refer to drawing 1 8 ) is treated, the B-spline 
functions B4 and j (x) are written as "Bj (x). ,! 

[0020] Then, spline-function F (x) is expressed like the following formula (3) as linear combination of 

the B-spline function Bj (x). 

[0021] 

[Equation 21 
n-f 

F(x)=:S PjB:(x) •"(») 

[0022] By asking for the coefficient pi of an upper formula (3), the spline function for which it asks 
becomes settled. 

[0023] Then, the simultaneous equations about a coefficient pi are stood. First, in each knot, if the value 

of input data and the value of a spline function are equal, it will set (refer to the following formula (4)). 

[0024] 

[Equation 3] 
n-1 

F(Xj) = .^: P i B i (x j )=y j , j=0, n -(4) 

[0025] To the number of Unknowns pi being n+3, since the number of input data is n+1, it runs short of 
two conditions by the upper formula (4). Then, the following formula (5) which are the endpoint 
conditions that the derivative of the second floor in ends is 0, and (6) are added. 
[0026] 
[Equation 4] 
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d 2 F 



d x' 



n-l 

x=x 0 



d 2 B 



d x 



= 0 



(5) 



x= 



x 0 



d Z F 



d x' 



n-l d z B j 
x=x n »" 3 



'dx^ 



= 0 



(6) 



x=x r 



[0027] pi can be found if an upper formula (4), (5), and (6) are solved as simultaneous equations. 
[0028] Also about the case more than two-dimensional, it can come back and solve to one dimension. 
Since the same is more said of high order origin, a two-dimensional case is explained. 
[0029] The inputted data shall be located in a line in the shape of a grid. That is, it considers as xi, i= 0, - 
-, m and yj, j= 0, --, the thing by which the value zij is given to (xi, yj) to n (refer to drawing 19 ). 
[0030] If extension of a knot is performed like [ direction / each / x directions and / of y ] the case of one 
dimension and the 3rd spline function is expressed using B-spline function, it will be given by the 
following formula (7). 
[0031] 

[Equation 5] ^ 

F (x. y) =jZj ^PjjBj (x) Cj (y) -(7) 

[0032] However, Cj (y) is B-spline function corresponding to yj. 

[0033] At this time, interpolation in one dimension is first performed to each i. That is, from input data 
zij, the function Fi of Variable y (y) is carried out like [ in the above-mentioned one dimension ], and it 
asks for it (refer to drawing 20 ). 
[0034] 
[Equation 6] 

n-l 

F. fv) —^la-.-.C: (y) —<8) 



i 



[0035] It performs 1 -dimensional interpolation each about j shortly to the solution qij of an upper 
formula (8). That is, it asks for the function Fj of Variable x (x) by the method in one dimension from qij 
(refer to drawing 21 ). 
[0036] 
[Equation 7] 



ij B i 



(x) 



[0037] pij obtained by the upper formula (9) serves as a coefficient of B- spline-function expression of 
the two-dimensional spline function for which it asks. 



[Translation done.] 



http://www4.ipdl.jpo.go.jp/cgi-bin/tran_web_cgi_ejje 9/4/2003 



Page 1 of 2 



* NOTICES * 

Japan Patent Office is not responsible for any 
damages caused by the use of this translation. 

1. This document has been translated by computer. So the translation may not reflect the original 
precisely. 

2. **** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[Drawing 1] It is drawing showing the composition of the 1st operation gestalt of this invention. 
[Drawing 2] It is a flow chart explaining operation of the operation gestalt shown in drawing 1 . 
[Drawing 3] It is drawing showing the composition of the 2nd operation gestalt of this invention. 
[Drawing 4] It is drawing showing the composition of the 3rd operation gestalt of this invention. 
[Drawing 5] It is a flow chart explaining operation of the operation gestalt shown in drawing 4 . 
[Drawing 6] It is drawing showing the composition of the 4th operation gestalt of this invention. 
[Drawing 7] It is drawing for explaining a principle and an operation of this invention. 
[Drawin g 8] It is drawing for explaining a principle and an operation of this invention. 
[Drawing 9] It is drawing for explaining a principle and an operation of this invention. 
[Drawing 10] It is drawing for explaining a principle and an operation of this invention. 
[Drawing 11] It is drawing for explaining a principle and an operation of this invention. 
[Drawing 12] It is drawing for explaining a principle and an operation of this invention. 
[Drawing 13] It is drawing for explaining a Prior art. 
[Drawing 14] It is drawing for explaining a Prior art. 
[Drawing 15] It is drawing for explaining a Prior art. 
[Drawing 16] It is drawing for explaining a Prior art. 
[Drawing 17] It is drawing for explaining a Prior art. 
[Drawing 18] It is drawing for explaining a Prior art. 
[Drawing 19] It is drawing for explaining a Prior art. 
[Drawing 20] It is drawing for explaining a Prior art. 
[Drawing 21] It is drawing for explaining a Prior art. 
[Drawing 22] It is drawing for explaining the trouble of a Prior art. 
[Drawing 23] It is drawing for explaining the trouble of a Prior art. 
[Description of Notations] 

101 Over Sampling Technique Means 

102 Filtering Means 

103 Interpolation Curvilinear Generation Means 

104 Error Calculation Means 

1 05 Input Means 

106 Output Means 

301 Over Sampling Technique Means 

302 Filtering Means 

303 Interpolation Curvilinear Generation Means 

304 Error Calculation Means 

305 Input Means 

306 Output Means 

401 Over Sampling Technique Means 
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402 Filtering Means 

403 Interpolation Curved- Surface Generation Means 

404 Error Calculation Means ( 

405 Input Means 

406 Output Means 

601 Over Sampling Technique Means 

602 Filtering Means 

603 Interpolation Curved-Surface Generation Means 

604 Error Calculation Means 

605 Input Means 

606 Output Means 

1301 Inputted Numeric Data 

1302 Point of Wanting to Calculate Corresponding Value 
1401 Partition Primary Function 

1501 Polynomial 
1601 Spline Function 
1801 3rd B-Spline Function 
2201 Spline Function 
2301 Spline Function 



[Translation done.] 
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[Drawing 7] 
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[Drawing 8] 




[Drawing 91 



[Drawing 10] 



[Drawing 2] 
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[Drawing 11] 



[Drawing 12] 
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[Drawing 51 
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[Drawing 22] 
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[Drawing 2\] 
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[Drawing 181 
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